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We describe current experimental and theoretical under- 
standing of the pressure-induced volume collapse transitions 
occurring in the early trivalent rare earth metals. General 
features of orbitally realistic mean-field based theories used 
to calculate these transitions are discussed. Potential defi- 
ciencies of these methods are assessed by comparing mean 
field and exact Quantum Monte Carlo solutions for the one- 
band Hubbard and two-band periodic Anderson lattice mod- 
els. Relevant parameter regimes for these models are deter- 
mined from local density constrained occupation calculations. 



I. INTRODUCTION 

Many of the trivalent rare earth metals undergo a dra- 
matic transformation in physical properties under com- 
pression, which is generally believed to arise from a 
change in the degree of 4/ electron correlation In 
some cases these changes appear abruptly across first- 
order phase transitions accompanied by unusually large 
volume collapses of 9 to 15%. Similar behavior is ob- 
served in the actinides, both for individual members un- 
der pressure, as well as across the series as a whole at 
atmospheric pressure Loosely speaking, the / 

electrons act as if they participate in the crystal bonding 
in the compressed, more weakly correlated regime, and as 
if they do not at larger volumes where correlation effects 
are more dominant. The terms itinerant and localized, 
respectively, are commonly used to describe the differing 
/ electron behavior in these two regimes. The intuitive 
concepts underlying these terms permeate the extensive 
investigations of both series of /-electron metals . 

Local density functional theory and its gradient ap- 
proximation improvements appear to do well for the more 
weakly correlated phases at pressures above these transi- 
tions, as may be judged by the considerable success ob- 
tained for Ce and the light actinides |^,^. Comparable 
predictive capabilities are lacking in the more strongly 
correlated regime, as is a satisfactory treatment of the 
large-volume collapse phase transitions themselves. 

It is the purpose of this paper to review current ex- 
perimental and theoretical understanding of these tran- 
sitions, with primary focus on efforts to develop a more 
rigorous and predictive treatment of strongly-correlated 
f-electron metals. One approach is the use of various 
corrected forms of local density functional theory which 
exhibit all-orbital realism, however, are still basically of a 
mean-field nature. General characteristics of these meth- 



ods are described and compared to Hartree Fock. In or- 
der to assess correlation effects neglected by these meth- 
ods, we report exact Quantum Monte Carlo (QMC) cal- 
culations for few-band effective Hamiltonians, and com- 
pare these to Hartree Fock results. In particular, new re- 
sults are reported for the three-dimensional two-band pe- 
riodic Anderson Hamiltonian, which represent a step to- 
wards using tera-scale computing resources to push QMC 
calculations into increasingly more realistic regimes. 

In the remainder of this paper. Sec. 2 reviews the 
relevant experimental data as well as some of the sim- 
pler theoretical concepts suggested by this data. Section 
3 provides parameter values which help place the rare 
earth transitions in the context of many-body effective 
Hamiltonians. Approximate but orbitally realistic the- 
oretical approaches to understanding the transitions are 
discussed in Sec. 4, while results of exact Quantum Monte 
Carlo solutions for few-band effective Hamiltonians are 
presented and discussed in Sec. 5. Our summary is given 
in Sec. 6. 



II. BACKGROUND 

This section reviews experimental data for the rare 
earth metals which characterize the issues of interest 
to this paper, as well as some of the associated calcu- 
lations. Important results include the pressure-volume 
equations of state, sequences of structural phase transi- 
tions, variations in equilibrium volume across the series, 
and the magnetic moments. The primary focus here will 
be on the trivalent rare earths in the first half of the se- 
ries, Ce (/i), Pr {p), Nd {p), Pm Sm and 

Gd {f), excluding divalent Eu {f). In the course of 
this review, it is useful to simultaneously discuss simple 
one-electron concepts which help to organize these re- 
sults. When viewed in this manner, one may interpret 
the data as showing evidence of /-electron participation 
in the crystal bonding in the higher-pressure more weakly 
correlated phases, whereas these manifestations appear 
largely absent in the lower-pressure more strongly corre- 
lated phases. 



A. Structural 

Figure |l| shows the room temperature pressure- volume 
curves for Ce [|, Pr g- Nd |l2|,|l3|, Pm Sm 
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FIG. 1. Pressure volume data for the rare earths. Struc- 
tures are identified, with "cmplx" signifiying a number of com- 
plex, low-symmetry structures. The volume collapse transi- 
tions are marked by the wide hatched lines for Ce, Pr, and 
Gd, while lines perpendicular to the curves denote the d-fcc to 
hP3 symmetry change in Nd and Sm. The curves are guides 
to the eye. Note that the data and curves have been shifted 
in volume by the numbers (in A^/atom) shown at the bottom 
of the figure. 

Ip^ - |]l7|| , and Gd , where some of the more recent, 

higher-pressure data is plotted, with curves provided as 
guides to the eye. More complete references may be ob- 
tained elsewhere from recent papers dealing with system- 
atics of the whole series For clarity of presenta- 

tion, the curves in Fig. |l| have been shifted in volume as 
shown at the bottom of the figure. The symbols change 
from open to filled or vice versa at structural phase tran- 
sitions, which as may be seen are quite frequent. The 
volume changes at these phase transitions are generally 
a few percent or less. There are three notable exceptions 
which are marked with the thick hatched lines at 0.9 
GPa in Ce (15% volume change) §, 20 GPa in Pr (9%) 
@~ 0, and 59 GPa in Gd (11%) These "volume- 
collapse" transitions demark the itinerant and localized 
regimes at pressures above and below the transitions, re- 
spectively, in these materials. Before discussing possible 
counterparts for the other three rare earths in Fig. |l|, it 
is useful to first review common structural trends which 
will be identified with the localized regime. 

The regular rare earths, those excepting Ce, Eu, and 
Yb, show a common structural sequence: hexagonal close 
packed (hep) — > Sm-type double-hexagonal closed 
packed (dhcp) — > face-centered cubic (fee) — > distorted 
fee (d-fcc) . The first four phases are all close packed and 
represent stacking variants of hexagonal layers such as 
the ABC ■ ■ ■ versus AB ■ ■ ■ order in fee and hep, respec- 
tively. There has been considerable discussion about the 
last phase |^^, which was initially named distorted- 
fcc due to the appearance of superstructure reflections in 
the fee diffraction pattern. The most recent opinion is 



that this phase is trigonal with eight atoms in the rhom- 
bohedral cell, and is due to a softening (TA) phonon at 
the L point in the fee Brillouin zone All or 

the latter part of this general hep —^ Sm-type —^ dhcp 
fee d-fcc sequence is observed under pressure in 
each of the regular rare earths. The valence electrons 
in these metals may be viewed as being compressed ei- 
ther by the application of pressure or by reducing the 
atomic number at fixed pressure p6| . Thus the lighter 
members of the series enter into the generalized sequence 
at successively later points, given their ambient, one at- 
mosphere, room temperature phases: hep (Gd), Sm-type 
(Sm), dhcp (Pm, Nd, Pr, and Ce). Aside from the low- 
est hep and dhcp regions for Gd and Ce, respectively, 
which lie below the range plotted, the remaining phases 
in the general sequence may be seen in Fig. Q for each of 
Pr-Gd, culminating in the high-pressure d-fcc end (filled 
up-triangles) . While both Pr and Gd are seen to undergo 
the volume-collapse transition at room temperature from 
the d-fcc phase, recent work has shown this phase to dis- 
appear above a 573-K triple point in Pr, so that the Pr 
collapse above this temperature occurs directly from the 
fee phase This behavior is similar to Ce, as seen in 
Fig. |l|, which undergoes the collapse transition at room 
temperature from the fee phase. 

The significance of the regular rare earth sequence lies 
in the fact that it appears to have no connection whatso- 
ever with / electrons. An elegant demonstration of this 
fact is provided by the occurrence of this structural series, 
including the d-fcc phase pO | , in compressed Y which has 
no nearby / states at all 27 . Theoretical calculations 



have furthermore demonstrated that the origin of the se- 
ries lies in increased occupation of the 5d states caused 
by pressure-induced shift in the relative position of 6sp 
vs 5d bands [||- ||]. 

The closed-packed and relatively high symmetry struc- 
tures of the regular rare earth sequence stand in sharp 
contrast to the low symmetry structures seen at pres- 
sures just above the volume-collapse transitions in Pr (or- 
thorhombic a-U structure) and Gd (body-centered mon- 
oclinic), both represented in Fig. m by open diamonds. 
While the Ce volume-collapse is isostructural, fee — * fee, 
at higher pressures in Fig. ^, one obtains a monoclinic 
or possibly the a-U structure (open diamonds), which is 
a subject of current debate Beyond this, Ce trans- 
forms at 12 GPa into the body-centered tetragonal (bet) 
phase shown by the filled down-triangles in Fig. |^, a phase 
also assumed by Sm above 91 GPa and similarly denoted. 

Low symmetry orthorhombic and monoclinic phases 
are of course prevalent among the early actinides, and 
have been demonstrated to arise from 5/ electron partic- 
ipation in the bonding P,^. Such phases are likewise 
taken to be evidence of itinerant 4/ character in the high- 
pressure rare earths. However, it should be noted that 
for low / band filling as well as for increasing pressure 
it is also possible to obtain higher-symmetry structures 
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in the presence of / electron bonding. The fee struc- 
tures of the collapsed a-Ce phase as well as Th [ ^3|j34| 
are characteristic of the metals, while high pressure 
bet phases are seen in Ce, Sm, Th, and predicted for U 
. Ultimately, U is predicted to reach an even higher 
symmetry bee structure, as are a number of the other 
early actinides Q . 

Finally, we turn to the three rare earths, Nd, Pm, and 
Sm, in Fig. |l| which do not appear to exhibit volume- 
collapse transitions. Recent experiments have identified 
Sm(V), indicated in Fig. |l] by the open diamonds, as a 
hexagonal phase with three atoms in the primitive cell 
(hP3). This is a rather unusual structure with open 
channels of roughly helical shape running along the c di- 
rection, and an approximate fourfold coordination |p|,|l7|. 
An analysis of anomalies in the pressure- volume equation 
of state for Sm has led to the conclusion that a rapid in- 
crease of 4/ bonding most likely begins with this phase 
[0. Thus in the last part of the Sm sequence, d-fcc 

hP3 bet, the 91 GPa transition between d-fcc and 
hP3, marked by the line perpendicular to the P-V curve 
in Fig. |l|, is tentatively identified as the onset of itinerant 
4/ character in Sm. The same d-fcc hP3 transition is 
seen in Nd |0, and is similarly marked in Fig. |[ Earlier 
work had noted a structural change in Nd at about 40 
GPa, denoted by the open diamond in the figure ||l2|,|l3| . 
More recent measurements report the hP3 phase from 
40 to 65 GPa in Nd, however do not publish pressure- 
volume data Data taken to date for Pm still falls 
within the regular rare earth sequence, as may be seen 
by the highest pressure d-fcc points in Fig. (filled up- 
triangles), so it is not yet possible to speculate as to the 
onset of itinerant behavior in this material. 



B. Equilibrium volumes 

Comparison of the atmospheric pressure or equilibrium 
volumes, Vq, of the / and d electron metals provides 
graphic evidence of regimes where the respective narrow- 
band electrons appear to be participating in the bonding 
or not, as has been appreciated for some years |^,^. 
Figure H shows these equilibrium volumes for the rare 
earths (4/), actinides (5/), and the 4d transition metal 
series [|7|. The horizontal axis is band filling, n/N, i.e., 
the nominal number of / or d electrons divided by the 
shell capacity, = 14 or 10, respectively. 

Aside from the two divalent rare earths, the depen- 
dence of Vq (n) for the rare earths is flat with a downward 
slope for increasing n associated with the lanthanide con- 
traction, i.e., the added / electrons do not fully screen 
the valence electrons from the increased core charge. In 
contrast, the Ad transition series has a parabolic shape, 
which may be understood from simple tight-binding ar- 
guments 1^ . Suppose the narrow band broadens symmet- 
rically with compression about the atomic energy level, 
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FIG. 2. Atmospheric pressure volumes of the 4d transition 
metals, the rare earths (4/), and the actinides (5/). The two 
rare earths with large volumes are divalent. 

and that it is characterized by a rectangular density of 
states, D{e), extending from —W/2 to W/2, as sketched 
in Fig. ^. The reduction in total energy, AE, associated 
with having a finite band width and the corresponding 
pressure correction, AP, is then 

AE^J deeD{e) = ---{N-n), (1) 

dE dW n 

AP = = (N-n), (2 

dV dV 2N^ ^ ' 

where is the Fermi energy. As the band width, W ^ 
grows with decreasing volume, V , Eq.(|^) implies a neg- 
ative contribution to the pressure which has a parabolic 
dependence on the band filling, n{N — n). Such a depen- 
dence is clearly in evidence for the 'id transition metal 
series. 

The actinides show mixed character as can be seen 
from Fig. ^. The first part of the series (Ac-Pu) ex- 
hibits similar parabolic dependence on band filling with 
increasingly lower values of Vq suggestive of a 5/ bonding 
contribution via Eq.(H). These same phases as already 
noted tend to exhibit complex low-symmetry structures 
also attributed to 5/ bonding effects. The latter part of 
the series (Am and beyond), on the other hand, acts like 
a second rare earth series in which the bonding contribu- 
tion, Eq.(0), no longer appears in evidence, and crystal 
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FIG. 3. Schematic rectangular density of states of a band 
of width W and capacity N , with /x the Fermi level. 
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structures are again of high symmetry . 

While the present paper is concerned with the rare 
earth transitions, the pronounced similarities between 
the two series of /-electron metals serves to strengthen 
the combined understanding. In both cases the valence 
electrons may be effectively compressed either by the 
application of external pressure or by the reduction of 
atomic number at fixed pressure [ p6| . Thus the drop 
in volume from Am to Pu at ambient pressure is the 
volume-collapse in the actinides, which may therefore be 
viewed as occurring at negative pressures for the first five 
members (Th-Pu) of the series. The rare earth series is 
simply off-set in pressure, requiring the application of 
positive pressure to drive any of the rare earths through 
the collapse. 

C. Magnetic moments 

Another key distinction between the localized and itin- 
erant regimes appears to be the presence or absence, re- 
spectively, of magnetic moments on the /-electron sites. 
This is an especially important signature for theoreti- 
cal simulations, since moment formation or loss is more 
amenable to study by simpler model calculations than is, 
for example, change in crystalline symmetry. Note that 
the loss of moment as detected by, e.g., magnetic suscep- 
tibility measurements may be due either to quenching of 
the /-shell moment itself or to its screening by the sur- 
rounding valence electrons. These two possibilities un- 
derly the Mott transition model of Johansson |3^, and 
the Kondo volume collapse model of Allen and Martin 
p9[ and Lavagna et. al. [Q, for the rare earth collapse 
transitions, respectively. In both cases it is important to 
understand when the /-shell can support a stable mo- 
ment. Simple one-electron concepts provide some insight 
into this issue, and are briefly discussed here after sum- 
marizing the experimental data. 

The localized phases with magnetic moments generally 
undergo some kind of magnetic ordering transition at low 
temperatures, with Gd (/'') having the highest ordering 
temperature at 293 K, followed by Tb (f) at 230 K 
among the rare earths pd|| , with the late actinides order- 
ing at ^^50 K or below |42| . Since the room-temperature 
data in Figs. ^ and |^ is all in the respective paramag- 
netic regimes, the critical characteristic is the moment 
itself and not possible ordering, which we shall ignore in 
this paper. As well be seen later in discussing the exact 
Quantum Monte Carlo calculations, magnetic order can 
in fact be an unwanted distraction. 

Figure ^ gives the observed atmospheric-pressure mo- 
ments in the rare earth and actinide ^] metals. The 
paramagnetic moments are extracted from the slope of 
the inverse magnetic susceptibility versus temperature in 
the paramagnetic regime, and may be compared to the 
corresponding free-ion values (solid line) given by 
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FIG. 4. Atmospheric pressure magnetic moments of the 
rare earth and actinide metals (symbols), and simple free ion 
estimates (curves). Both paramagnetic and saturation mo- 
ments are shown. 

g[J{J + In the latter, one assumes Russel- 

Saunders coupling and takes the Hund's rules ground 
state for each ion to obtain the total angular momen- 
tum, J, while g is the Lande factor. Saturation moments 
are obtained either from the value of the magnetization 
at low temperature, or from extrapolations of magneti- 
zation data to infinite field. The corresponding free-ion 
value here is gJ as plotted by the dashed line in Fig. ^. 

It is evident from Fig. ^ that atoms in the late 
rare earth and actinide metals exhibit moments very 
close to their free ion values, which attests to / elec- 
trons relatively unperturbed by their crystalline envi- 
ronment, consistent with localized / electron behavior. 
The early actinides, Th-Pu, on the other hand, all ex- 
hibit temperature-independent paramagnetic susceptibil- 
ity consistent with the absence of magnetic mo- 
ments. This is a dramatic change from free-ion be- 
havior suggesting strong interaction between the / elec- 
trons in these metals and their local environment. A 
temperature-independent paramagnetic susceptibility is 
also observed for Am |Q, however, the /^, J — ion 
should not have a moment even in the free ion limit, so 
that magnetic behavior is not particularly illuminating 
in this case. 

The early rare earths are somewhat more complicated, 
since crystal field interactions split their Hund's rules 
multiplets impacting their magnetic properties. The sat- 
uration moment of 2.7 /ie for Pr in Fig. ^was obtained 
at 4.2 K for example, while at a temperature on the 
order of the crystal field splitting, 40 mK, the size of the 
moment obtained from neutron scattering appears to be 
0.36 /iB [§■ Nevertheless, these crystal field splittings 
are still rather small, so that the 2 J-fl-level multiplet and 
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its associated magnetic moment may be viewed as intact 
on the more coarse energy scale of interest to this paper. 
These crystal field interactions do serve as a reminder, 
however, of an approaching threshold beyond which the 
crystalline environment will destroy the local moments. 

Simple one-electron concepts provide some under- 
standing of the competing factors which determine 
whether an /-shell can support a stable moment. The 
Coulomb interaction, /, within the manifold of / states 
on a particular atomic site, j, may be written 

/ ~ fiif (fiif — 1) + multipole terms , (3) 

where hij = X^ma ^Ifma^ifma IS the total number opera- 
tor for all 14 / states indexed by ma, and the monopole 
Slater integral, = [//, is the usual scalar Hubbard pa- 
rameter describing the repulsion between each pair of / 
electrons on the site. A Hartree-Fock expectation of the 
monopole part of Eq. yields 

{I)o = \F%{trp^fY-trpif], (4) 

= + (5) 

where riif = tr pif, and the riia are the eigenvalues of the 
// one-particle density matrix for the site i, pi/, given by 

{Pif)ma,m'(j' = {^ifma^if™-''^') ' (^) 

It is the second term in Eq.(^) which is interesting, as 
it is minimized by having "sharp" orbital occupations, 
TT-ia = or 1. This is achieved in Hartree-Fock by one- 
electron energies, eka, containing the term Uf{nf —Ua), 
where n/ and Ua are site- independent values of n^/ and 
riia , respectively. This term discriminates between occu- 
pied and empty states, placing the latter higher by Uf 
when the Ha = 0, 1. The multipole terms in Eq.(|^) sup- 
plement Eq.(||) in impacting the nature of the orbitals 
which diagonalize the // density matrix, p, as well as 
helping to select those which will be populated. The 
exchange terms, for example, favor all parallel spins for 
n-if < 7 (Hund's first rule). Other multipole terms seek 
to maximize the total angular momentum for the / shell 
(second rule), or an approximation to this in mean field. 
The net effect in the exact (with intraatomic correlations) 
solution is of course to build large values of the total spin, 
Sif, and total angular momentum, Lif, from the multi-/- 
electron shell, which may then be combined by spin-orbit 
(third rule) into a total moment gJ. 

This behavior is in contrast with the one-electron 
bonding effect in Eq.(0), where to a rough approximation 
one can find all spin and orbital types in the lower part 
of the band. This follows from the fact that the overall 
size of the one-body, intersite hopping matrix elements, 
{4'ifim\H\(f>jfim'), is set by the tail of the radial wave- 
function which to first approximation is independent of 



m, and that geometric effects on these matrix elements 
tend to be smeared out both by the variety of neighboring 
atom positions as well as by the range of phase relations 
associated with points throughout the Brillouin zone. As 
a consequence, the band broadening mechanism of Eq. (|^) 
is relatively indiscriminate in terms of spin-orbital occu- 
pations, favoring roughly n.if/lA which serves to 
quench spin and orbital moments. 

It seems clear in the extremes Wf » Uf and Wf « 
Uf that the competition between Eqs.(^ and (||) results 
in the absence or presence of a moment, respectively. It is 
also evident how at least Hartree-Fock goes about reduc- 
ing the / bonding effect of Eq.([|) in the latter regime, by 
splitting the occupied and empty states to create a new 
entirely full band for which n{N — n) — 0. The major 
current debate concerns how the experimentally observed 
moment is actually first lost with increasing pressure (or 
Wf/Uf) for the 4/ and 5/ metals, whether the / mo- 
ment is first quenched similar to the manner discussed 
here and by the Mott transition model |^^, or whether 
while still in a robust state it is first screened away by a 
surrounding cloud of valence electrons as in the Kondo 
volume collapse model p9| , ^ . 

III. PARAMETERS 

To provide contact with the intuitive concepts dis- 
cussed in the previous section as well as to motivate 
model problems which can be exactly solved by Quan- 
tum Monte Carlo techniques, it is useful to characterize 
the rare earth metals in terms of the effective Coulomb 
repulsion between / electrons on the same site, Uf, the 
position of the / level, ef, and the /-/ and /-valence 
hopping interactions, tff and tfy, or the equivalent ex- 
pressions in band widths, Wff and Wfy, respectively. 
These results were obtained as a function of compression 
for fee phases of the rare earths using a scalar-relativistic 
linear muffin-tin orbitals (LMTO) method in the atomic- 
sphere approximation plus combined correction |^ , ^ . 
All electrons were treated self consistently. The specific 
approximations used to calculate the different parameters 
have been discussed at length elsewhere . While 

not rigorous, they are reasonable, and known to give val- 
ues of the Coulomb interactions and hopping parameters 
within ^ 20% of values deduced from experiment in other 
materials ||,||,||. 

The monopole contribution to the /-/ Coulomb inter- 
action, Uf, and the site energy, Sf, may be obtained from 
self-consistent calculations of the total energy, E{n), or 
the /-orbital eigenvalue, A/(n) ~ dE{n) / dn, as a func- 
tion of / occupation, n. Although more sophisticated 
methods exist for decoupling the orbitals so as to control 
their occupations [Q , the 4/ orbitals are sufficiently well 
localized in the present materials that we have simply 
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treated them as part of the self-consistent core |Q . The 
formal definition of Uf is 

Uf^E{nf + l) + E{nf~l)~2E{nj), (7) 
«A/(n/ + i)-A/(n;-i), (8) 

while the site energy, Sf, may be defined in terms of the 
removal energy, E{nf) — E{nf — 1), 

ef + Ufinf-l)^E{nf)-E{nf-l) (9) 
^Xfiuf-^). (10) 

These equations provide the best quadratic approxima- 
tion, Ea + EfU Ufn{n — l)/2, to E{n) in the vicinity 
of the nominal integer occupation, n = nf. In principle 
one should probe the dependence on occupation of only 
a single / site in the infinite solid. For metals, however, 
the screening is so effective that one may in practice al- 
ter the occupation at every site. This requires use of a 
neutralizing electrostatic background, or exchanging the 
electrons between the / site and the valence band Fermi 
level, /i. In the latter case, case — /i should be added to 
the left side of Eq.(|). 

Figure |^ shows Uf versus atomic volume for the five 
rare earth metals which have exhibited anomalies un- 
der pressure. The pluses mark the boundaries of the 
volume collapse transitions (Ce, Pr, Gd) or the d-fcc 
hP3 symmetry change (Nd, Sm). The open squares 
are the atmospheric-pressure values calculated by Herbst 
and Wilkins however, omitting their correlation and 
Hund's rules corrections to be consistent with the present 
work. These corrections to Uf are less than 1 eV for Ce, 
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FIG. 5. Calculated Coulomb interaction, Uf, for the rare 
earth metals as a function of volume. 
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FIG. 6. Calculated removal energy, £/ -I- f//(n/— 1), relative 
to the valence Fermi energy, for the rare earth metals as 
a function of volume, n/ is the nominal integer / occupation. 

Pr, Nd, and Sm, however, increase the Gd value by 5.4 
eV. With the corrections, Herbst and Wilkins found quite 
good agreement with electron spectroscopy data [ p9[ . 
The general trends evident in Fig. || are as expected. The 
heavier rare earths have more narrow 4/ bands, reflect- 
ing more compact orbitals, and therefore exhibit larger 
values of Uf. In all cases the effect of compression is to 
enhance screening and therefore reduce U f . In the range 
of interest these reductions are ~ 20%. 

The removal energies, E(nf ) — E(nf—\) k, e f + Uf{nf— 
1), are shown in Fig. ^ relative to the Fermi energy, fj,y, 
of the three valence electrons per site. As before, the 
pluses mark the location of the anomalies, while the open 
circles simply locate the atmospheric-pressure volumes. 
The comparable values of Herbst and Wilkins (ab- 
sent correlation and Hund's rules corrections) show about 
the same dependence on atomic number, however, are 
^ 2 eV higher in energy than the present results. Their 
correlation and Hund corrections then lower these values 
by 1-3 eV for Ce-Sm (4.6 eV for Gd), so that their final 
results are within about an eV of the uncorrected values 
in Fig. H for Ce-Sm. Herbst and Wilkins found the same 
characteristic volume dependence as also seen here. In 
particular, the maxima seen in Fig. ^ in the range 15- 
20 A'^/atom corresponds to the end of an electronic s-d 
transition caused by the rising 6sp levels dumping their 
electrons into the 5d states. As they note, the shifting 
balance of kinetic and potential energies under compres- 
sion will eventually cause the 5d band to move upward 
relative to 4/ states, as seen at the lower volumes. The 
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FIG. 7. Calculated ratio, t///W^//, where W;j is the / band 
width due to /-/ hybridization, for the rare earth metals as 
a function of volume. These results were obtained for an fee 
lattice, so that near-neighbor hopping interactions, may be 
estimated by dividing the corresponding band widths by 16. 

important conclusion to draw from Fig. |, is that the 4/ 
levels are rising as a function of compression throughout 
the range leading up to the pressure anomalies. 

The contribution to the 4/ band width due to // hop- 
ping interactions may be simply determined by the width 
of the pure / Wigner-Seitz band, which is given by the 
energy difference between / orbitals whose logarithmic 
derivatives at the atomic-sphere boundary are (band 
bottom) and infinite (band top), respectively. The 4/ 
band widths, W//, obtained in this way for the five rare 
earths are plotted in Fig. |^ in the ratio J7/ /H^// as a func- 
tion of volume. As before, circles and pluses mark values 
at the atmospheric-pressure and the transition volumes, 
respectively. The heavier rare earths may be expected 
to have both narrower / bands and larger Coulomb in- 
teractions, leading to the ordering shown. The functions 
[//, WJ^, and UfjWff depend on the 0.2, 2.0, and 2.2 
powers of volume in compression, to within about ±0.05 
in these exponents. Canonical band arguments suggest a 
d^'^' ~^ dependence for the overlap of / and V orbitals a 
distance, d, apart or ~ F^ '^, which is close 

to what is observed here. 

The width in energy over which / states are spread also 
arises from /-valence hopping interactions or hybridiza- 
tion. We have obtained root-mean-square widths, W^''™^, 
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FIG. 8. Calculated ratio, W]'^'' IW]^ 
metals as a function of volume. Here 
rms / band widths due to /-/ and /-valence hybridization 
respectively. 

using / partial densities of states, Df{£), obtained in 
three ways. In one case we have taken the orthogo- 
nalized one-electron Hamiltonian matrices obtained from 
the self-consistent one-electron potentials and set all /- 
valence matrix elements to zero, so that the resultant 
band width, Wj™^, is due predominantly to // inter- 
actions with some small crystal field contributions. In 
the second case, we replace the // blocks by the identity 
matrix times the average / energy, obtaining a width, 
WJ^^, due to /-valence interactions. In the third case, 
we have used the unmodified orthogonalized Hamiltonian 
matrices, obtaining the full / width, WJ"^^. Note that it 
is rigorously the case for these rms widths that 



11/2 



(11) 



(12) 

The ratios, Wj^^/Wjp^, are shown in Fig. ^ for the five 
rare earths. In compression they depend on the 0.16 
power of the volume to within ±0.06 in the power. If the 
/-/ overlap acts like I ~ 2.5 in its d"'"' ~^ dependence 
{Wff ^ T^~^ °), then the volume dependence of the ratio 
in Fig. ^ suggests a valence / ~ 2. This is consistent with 
the general belief that the important intersite /-valence 
hybridization is with d states. 

The implication of Fig. || is that /-valence hybridiza- 
tion is a more important contributor to the overall / band 
width in the rare earths than is direct /-/ hybridiza- 
tion. Earlier work suggested comparable impact, based 
on the boundaries of one-electron bands of / symmetry 
pot . Equation ( pdj ) , however, also measures / orbital ad- 
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mixtures which can occur farther out in energy as parts 
of bands identified by symmetry as of other I character. 

If we define Wfy = Wff Wf^yWJJ^'', then Figs. |-| 
imply that the rare earth collapse transitions (or sym- 
metry changes in Nd and Sm) occur over the ranges 
Uf/Wfy = 3.7-2.6 (Ce), 1.8-1.5 (Pr), 1.6 (Nd), 1.1 (Sm), 
and 1.9-1.5 (Gd). These ratios are all in the range 1-2 
except for Ce, which has the lowest band filling. If 
the combined / band width, W/, is given by an expres- 
sion such as Eq.(p^), then the corresponding values for 
Uf/Wf would be reduced by 5%-7%. 

To make contact between the parameters presented in 
this section and the models to be discussed in Sec. 5, note 
that for near- neighbor interactions and only tff hopping, 
the full / band width, Wj = Wff, is 16 1// for the fee 
structure used in Fig. while 12 i// for the simple cubic 
structure assumed in Sec. 5. As a matter of convenience, 
and to acknowledge the scaling in Eq.(p^), we shall also 
define the width Wf — Wfy = 12 tfy for the simple-cubic 
two-band periodic Anderson model calculations in Sec. 5 
with a dispersionless (t//=0) / band. 



IV. MODIFIED MEAN FIELD THEORIES 

It is generally believed that the volume collapse transi- 
tions also occur at zero temperature, in which case they 
should in principle be reproduced by density functional 
theory |5l] as properties of the ground state total energy. 
In practice, however, local density functional (LDF) the- 
ory, which resembles a modified Hartree mean field [ p2[ , 
does not give the transitions. There has been much effort 
to find improvements which are more successful in this re- 
gard. These include the use of spin |Q- |5^ and orbital 
1^ , ^ polarization, self-interaction corrections [|6ll , 
and the LDA-I-U method A critical comparison of 

these methods has been given recently for the related 
problem of transition metal monoxides []63| . Here, we 
briefly review such calculations for the / electron metals. 
Then we describe results of screened Hartree-Fock calcu- 
lations on orbitally realistic effective Hamiltonians which 
are similar to some of these methods, and by this means 
illustrate a number of common features of the modified 
mean field theories. 



A. Corrected local density functional theory 

One of the most successful applications of spin- 
polarized LDF theory was to the actinides, in which the 
jump in equilibrium volumes seen in Fig. |2| was repro- 
duced , as well as the delocalization transition in com- 
pressed Am p4|| . Similar calculations for compressed Ce 
p5[ and Pr p6| were less successful, with too little of 
the / bonding contribution removed in the large volume 
spin-polarized region. Note that the collapse between Pu 



and Am in Fig. || occurs near half filling of the 5/ shell 
so that spin polarization potentially reduces n[N — n) in 
Eq.(||) by a factor of eight for Am, in contrast to only 
about a factor of two for the early rare earths, Ce and 
Pr. 

The solution to this difhculty has been orbital polar- 
ization [ prf , in which a term proportional to the square of 
the z component of the total angular momentum is added 
into the total energy functional so as to simulate Hund's 
second rule. Combined with spin polarization, the result- 
ing calculations discriminate amongst the 14 / states in 
favor of large values of the z components of total spin, 
Ms, and angular momentum. Ml, for the / shell. Both 
effects arise from the multipole terms in Eq. (^ , and their 
size is dictated by combinations of the material specific 
values of Slater's F"^, F^, and F^ integrals. The conse- 
quence is the possibility of splitting off occupied bands 
of any integral number of / electrons per site, i.e., re- 
ducing n{N — n) to potentially zero for any value of n. 
Such calculations for the collapse transitions in Ce |5^] 
and Pr have yielded reasonably good agreement with 
experiment, in the latter case supplemented also by gen- 
eralized gradient corrections to the exchange-correlation 
potential. Plots of the spin and orbital polarization in 
the first case show roughly maximal values in the local- 
ized regime, which decay to zero in the vicinity of the 
transition | |57[ |. 

One may also discriminate amongst the 14 / orbitals 
purely on the basis of their occupation, as noted in regard 
to the second term in Eq.(||), with the size of the effect 
given by the monopole integral, F^ , the familiar Hubbard 
U f. Because of the poorly cancelled self-interaction of an 
electron with itself in LDF theory as well as its spin and 
orbitally polarized generalizations, this effect is largely 
absent from the above calculations. It may be rein- 
troduced by performing self-interaction corrected (SIC) 
LDF calculations generally combined with spin polariza- 
tion. Such calculations split the occupied and empty / 
states, determined self-consistently, by Uf and therefore 
also serve to reduce n{N — n) of Eq. (^) potentially to 
zero, for any value of n, in the large volume localized 
regime. Quite satisfactory results have been obtained for 
both the Ce and Pr volume collapse tran- 

sitions. The LDA-I-U method ||64| is an ad hoc way of 
achieving the same end in a much easier calculation by 
simply adding the second term in Eq.(^) to the LDF- 
total energy functional. Apparently, a value of J7/ 3 
eV, about half of what would be expected, is required 
to make the Ce transition occur in the right place [ |62[ . 
This is interesting, since judging from the separations be- 
tween the centers of gravity of the empty and occupied 
4/ states, SIC calculations for Ce H and Pr |6|] both 
correspond to Uf ~ 9-10 eV. It might also be noted that 
~40% smaller effective U values are also required in ana- 
lytic random phase and conserving approximations when 
attempting to fit exact Quantum Monte Carlo results 
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B. Hartree Fock with static screening 

The Hartree-Fock (HF) method has no self interac- 
tion problem. Unfortunately, its electron-electron inter- 
actions are unscreened so that the one-electron spectrum 
reflects the full value of the bare Slater monopole inte- 
grals, e.g., F" ~ 26 eV for the 4/ states in Ce 
If the Hamiltonian is written in second quantized form, 
however, it is trivial to replace F° by a screened value. 
HF solution of such a Hamiltonian includes some correla- 
tion effects via this static screening, and might be viewed 
as a further approximation to the GW method ||6^ . 

We report here such Hartree-Fock calculations for the 
Hamiltonian, 

k,Q,/3 k,m,(T 

+ ^UfY,n^f{n^f-l), (13) 

i 

where k are vectors in the Brillouin zone, a = Ima 
represents the usual angular, magnetic, and spin quan- 
tum numbers, i are the lattice sites, and fiif — 
12ma ^Ifma'^ifina IS the total / number operator for site 
i. We consider only one atom per unit cell, so that 
the localized HF solutions of Eq.(p^) will be ferromag- 
netic. Rather than the limited hopping parameter infor- 
mation given in Figs. 0-^, we use the full, converged, k- 
dependent LDF Hamiltonians, /i)^J^^(k), which are 32x32 
matrices {s-f basis plus spin) that have been orthogo- 
nalized. They are corrected by replacing the LDF / site 
energy by its constrained occupation counterpart from 
Fig. y, and adding the screened, /-/ monopole Coulomb 
interaction, though as yet, no multipole contributions 
such as exchange. Note that each of h^^^ , , £f, 
and Uf is volume and material dependent. 

Figure |^ summarizes the ground state HF solutions of 
Eq.(|l3|) for Ce. We have omitted the spin-orbit interac- 
tion here for simplicity, however, we obtain similar results 
with this interaction added to Eq. (|l^) . The figure shows 
the total / occupancy per site, Uf — riif = (fiif) as a func- 
tion of position of the / level, £/, relative to the valence 
Fermi level, ^v, i-e. that for just three spd electrons per 
site. Each curve corresponds to a different volume. The 
calculated values of e/ — /i^ from Fig. ^ are shown as the 
open circles on the appropriate curves. The use of e/ as 
an independent variable here is both instructive and of- 
fers some sense of the impact of the ~ 1 eV uncertainty 
in these values discussed in Sec. 3. 

Note first that at large volume (solid curve, V — 
56.7 A'^ /atom), where the hybridization interactions are 
weak, one obtains a staircase structure of integer values 
for Uf. As £/ drops below /i^, only one / electron state 




Ef-Hv (eV) 

FIG. 9. Static-screened Hartree-Fock calculations for Ce. 
The number of / electrons per site, n/ , is shown as a function 
of the / level position, Sf, relative to the valence Fermi energy, 
fly , at four different volumes. The calculated values of from 
Fig. 1^ are shown by the open circles. The solutions exhibit 
different branches characterized by integer numbers, n^"*^, of 
occupied split-off / bands. 

is occupied at first due to the cost, Uf, for adding a sec- 
ond. Only when Sf is lowered an additional Uf below /z„ 
is the second / electron picked up, and so on. Such a plot 
is similar to an integral over the /-decomposed density of 
states, so that the plateaus correspond to Mott-like gaps 
in the / spectrum. A related plot, where the chemical 
potential is varied, is a standard tool used in analyzing 
the results of many-body simulations |Q . This staircase 
structure is smeared out with the increasing hybridiza- 
tion width of the / band as volume is reduced. When 
this width becomes so large that Uf is effectively unim- 
portant, Uf grows linearly with decreasing ef — fj,y. The 
idealized limits of localized and itinerant character are 
associated with such staircase versus linear behavior, re- 
spectively. The usefulness of a plot such as Fig. ^ is that 
it provides a visual sense of where a given material at a 
specific volume lies between these two limits. In partic- 
ular, the short dashed curve in Fig. ^ corresponds to a 
volume very close to that of the collapsed a phase of Ce, 
suggesting that the Ce transition lies closer to the local- 
ized limit. The ratios Uf/Wff = 6.7 and Uf/Wfy = 2.6 
from Figs. 0-^for Ce at the a volume are consistent with 
this observation, and might be compared to 3-5 and 1.1- 
1.5, respectively, on the low- volume side of the other rare 
earth transitions. 

We have labeled the HF solutions in Fig. ^ by integer 
values of n'^'^ , referring to the number of / electrons split 
off in low energy fully occupied bands. The partial / 
density of states (DOS) shown in Fig. |l^ illustrates the 
differences at essentially the a-Ce conditions {V — 27.8 
A'^ /atom and e f—^J-v = —2.2 eV) . The ground state n^'^ = 
1 solution (dashed/shaded curve) exhibits the features 
associated with one localized electron: 98% occupancy of 
a single spin-orbital (total n/ = 1.05), and a separation 
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FIG. 10. Ce / densities of states for stable v}f"^ = 1 and 
metastable n^f"' — Hartree-Fock solutions for Ce at essen- 
tially the a phase volume and £/. The zero of energy is the 
full system Fermi energy, fi. 

between occupied and empty / states of ^ [// = 6.7 cV. 
The metastable n'^"^ = solution, on the other hand, has 
its rif = 0.71 / electrons distributed across all of 14 / 
states, and is consistent with itinerant / electron behav- 
ior. While one might have expected the latter solution to 
have the lower energy, Sandalov et. al. have found a 
similar result in their tight-binding and LDA-t-U calcu- 
lations, observing that they would need to reduce Uf to 
^2.7 eV to put the itinerant state lower. 

While it occurs at too high compression. Fig. ^ does 
ultimately predict a transition in Ce, since the boundary 
between n'^'^ = 1 and solutions shifts to lower values 
of Ef — fly with compression. The open circle for V = 
19.4 A'^/atom (dash-dot curve) may be seen to lie on the 
n^'^ = branch. We find the same behavior in compressed 
Pr however, the analogous transition is from n^ji"^ = 2 to 
1, as if the two / electrons in Pr delocalize one at a time 
in the present calculations. 

The present HP results could be improved by including 
multipole terms and, as noted by Sandalov et. al. [ p2[ , 
by using state-dependent hopping interactions. Nonethe- 
less, the present results illustrate some features which 
appear to be characteristic of the modified mean-field 
theories as a class. In all cases these theories model 
the localized regime by polarized solutions with an in- 
teger number, n^f"^, of filled split-off bands of relatively 
select spin-orbital character. The pressure induced tran- 
sitions predicted by these methods involve a shift of these 
bands to the vicinity of the Permi level, i.e., a transfer 
of spectral weight in roughly integer units of electrons 
per atom. Moreover, these transitions are abrupt in the 
sense that there are generally two distinct solutions (e.g.. 
Pig. |l^) with total energies which cross at some point. 
To anticipate the next section, correlation effects enable 
a more continuous and nonintegral transfer of / spectral 
weight from the low-lying localized or polarized / states 



to higher energies |^l|. In the case of Ce, however, Svane 
[ ]6l| has claimed that such effects modiiy the total en- 
ergy in the vicinity of the crossing, but not further away 
where the Maxwell common-tangent touches the two en- 
ergy curves. 

V. EXACT QUANTUM MONTE CARLO 
CALCULATIONS 

As emphasized above, electronic correlations play an 
important role in volume collapse transitions, yet ap- 
proximate mean-field treatments of their effects, such 
as Hartree-Pock theory, suffer serious deficiencies. We 
therefore turn to an approach, Quantum Monte Carlo 
(QMC), which can treat electron-electron interactions 
exactly. QMC is computationally very expensive, and 
while models of increasingly accurate orbital realism can 
now be simulated, it is still not feasible to study a sys- 
tem with the full complexity of, for example, Ce. A cru- 
cial task is to put together LDP calculations which do 
treat the full electronic structure, with a QMC treatment 
of simpler Hamiltonians which can focus on the physics 
missed by LDP theory. 

A. Quantum Monte Carlo method 

The idea of QMC is to use an exact mathematical 
transformation to rewrite the interactions between two 
electrons as a coupling of a single electron with a classical 
fluctuating field, which then mediates the interactions in- 
directly. The resulting single-particle quantum mechan- 
ical problem can be solved, leaving only an integration 
over all possible values of the classical field. This inte- 
gration is in an extremely high dimensional space. Por 
every pair of interacting orbitals in the original Hamil- 
tonian, there are L auxiliary field variables, where L is 
roughly the ratio of the bandwidth to the temperature. 
Such high dimensional integrals can be done by classical 
monte carlo methods. 

This formalism is straightforward and exact, but there 
are a number of difficulties in its implementation. The 
first is the scaling of the computation time with system 
size. While the integral to be done is classical, the inte- 
grand is non-local. Each field variable interacts with all 
the others, not just with a few "neighbors". To update 
every component of the field takes a computation time 
which scales as the cube of the number of orbitals. Thus 
to change the number of orbitals per site from a Hamil- 
tonian with one conduction and one local orbital per site, 
to a fully realistic all-orbital description in which 16 or- 
bitals are retained, increases the computational expense 
by a factor of 8"^. 

Going to lower temperatures also involves increased 
effort, since the number of field variables scales as band- 
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width over temperature. This is only a linear cost. That 
is, to double (3 = 1/T only doubles the cpu requirements. 
However, a more serious difficulty is also encountered. 
In many situations the integrand, which one is using as 
a sampling probability, can go negative. When this oc- 
curs, QMC simulations are no longer feasible. This "sign 
problem" does not occur for certain "symmetric" choices 
of the parameters in the Hamiltonian, and also does not 
occur at high temperatures. Nevertheless, it is a signif- 
icant limitation to the fillings and interaction strengths 
for QMC in general. 

The state-of-the-art in QMC on vector supercomput- 
ers are simulations of Hamiltonians with roughly a hun- 
dred orbitals. This means periodic lattices of, for ex- 
ample, 4x4, 6x6, 8x8, and 10x10 for a two dimensional 
single orbital model, and 4x4, 6x6, and 8x8 for a two 
dimensional two orbital model |f7^-[74|. Such sets allow 
for extrapolations to the thermodynamic limit, particu- 
larly when the form of the finite size correction is known 
|75| , [76| . It is clear that simulations with increased orbital 
realism, in three dimensions, require significantly more 
powerful hardware, such as is now becoming available 
with the parallel platforms. QMC is ideally suited for 
parallelization; it achieves an almost linear speedup as 
more and more processors are applied to a problem. We 
will discuss finite size effects in more detail later, but for 
the moment let us say that we do not expect them to be 
too significant for the thermodynamics, since the energy 
involves only local quantities [Q. Therefore our work 
will focus on the issue of more orbitals per site, as op- 
posed to increased spatial size. This is also the direction 
desired for increased contact with LDF theory. 

B. The Hubbard and Anderson Lattice Hamiltonians 

The determinant Quantum Monte Carlo calculations 
we will use consider electrons which move on a discrete 
lattice of atomic sites, and their associated orbitals WA - 
Continuum QMC methods, like Green's Function Monte 
Carlo, do exist, but as yet appear less easily applied 
to the problems of magnetic moment formation, mag- 
netic ordering, and singlet formation, which are relevant 
to the volume collapse transition. The simplest lattice 
model which might be applied to the present system is 
the single-band "Hubbard Hamiltonian." Originally for- 
mulated for "d" electron systems, in the present case one 
might consider a single effective / electron. If we denote 
by an operator which creates such an electron on site 
i with spin cr, then this Hamiltonian may be written 

H = -tff E (flfi- + 4/-) 

+Uf^{nii -^){nii-^) - H^n,^ . (14) 

i ia 



The kinetic energy and interactions in this model are as 
local as possible, hopping between near-neighbor sites 
{ij) only, and interactions between the density of elec- 
trons, rii^ — fj^fia, only on the same site. The interac- 
tion term has been written in a special form which makes 
"half-filling," a density of one electron per site on aver- 
age, occur precisely at fi = for any choice of temper- 
ature T or parameters tff,Uf in the Hamiltonian. Fur- 
thermore there is no sign problem. These desirable prop- 
erties require a bipartite lattice, composed of two sublat- 
tices, such that all near neighbors of an atom on one sub- 
lattice lie on the other. Body-centered cubic and simple 
cubic are examples, however, face-centered cubic is not. 
Since, for near-neighbor interactions, the first structure 
has a divergent noninteracting density of states at fi = 0, 
the present simulations have considered the simple cubic 
structure. 

The kinetic energy term in Eq. ( |l^ can be equivalently 
written in terms of operators fl^ which create electrons 
of momentum k as '^^^ ^kafl^fka- For a 3-d simple 
cubic lattice 

efco- — ^2t f f (cos k^a + cos kyU + cos k^a) , (15) 

where a is the lattice constant. This gives a bandwidth 
Wf = \2tff. As a function of atomic volume, V, Figs. ^ 
and shows the parameters in Eq. ( p^ to scale roughly 
like Uf ~ constant and t// oc Wf ^ V^^. We shall there- 
fore present the phase diagram of the one band Hubbard 
model in dimensionless units T/Uf versus Uf/Wf to ap- 
proximate T versus V^, respectively. 

Despite its simplicity, even this single band Hubbard 
model cannot be solved analytically, and while some of 
its properties are fairly well established, other features, 
both qualitative and quantitative, are subject to consid- 
erable debate j?^. What is fairly certain is that on a 
3-d simple cubic lattice the ground state is an antifer- 
romagnetic insulator at half-filling, for all values of the 
ratio Uf/Wf. The magnetic ordering can most simply 
be understood as a consequence of the divergence of the 
noninteracting magnetic susceptibility, 

. 1 tanh(/?£p/2) 
Xo(q*,r) = -^ ' (16) 

at ordering wavevector — (tt, tt, tt). Since Xo(q*j T) 
oo as T — > 0, the Stoner criterion UfXa = 1 is satisfied for 
any value of the interaction strength at sufficiently low 
temperatures. At strong couphng, large Uf/Wf, the in- 
sulating behavior is usually explained by noting that for 
densities less than half-ffiling the system can be in con- 
figurations in which no sites are doubly occupied, that is, 
no sites contain both an up and down spin electron. How- 
ever, at densities greater than half-filling, double occu- 
pation is inevitable, and hence the system suddenly finds 
itself in a manifold of states of an energy Uf higher. Thus 
the density of states of the model consists of "upper-" 
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FIG. 11. Phase diagram of the single-band Hubbard 
Model; 1 effective f orbital per site. The Neel tempera- 
ture Tjv below which antiferromagnetic ordering takes place 
is indicated by the boundary between the shaded (antiferro- 
magnetic) and unshaded regions. For small Uf/Wf, Tn is 
well-approximated by Hartree-Fock (solid line). For large 
Uf/Wf, Tn matches onto the spin-1/2 antiferromagnetic 
Heisenberg Model (dashed line). QMC, which is exact, picks 
up both limits and the intermediate coupling regime | pot , 
which both Hartree Fock and Heisenberg fail to capture. 
The inserts show the density of states for the infinite dimen- 
sional Hubbard model, and the continuous metal-insulator 
(MI) transition where the central resonance is lost |Q. The 
solid curve below is an extension of this MI boundary which 
might apply were the magnetic order suppressed 

and "lower-Hubbard bands" , which are separated in en- 
ergy by the on-site repulsion Uf. If C// exceeds the band- 
width of the system, this gives rise to a "Mott-Hubbard" 
gap at half-filling. At weak couphng, small Uf/Wf, these 
upper and lower Hubbard bands overlap, but the system 
remains insulating because of the doubling of the unit 
cell due to the antiferromagnetic order. 

Considering now finite temperature, but remaining at 
half-filling, one finds that the Neel temperature T/v be- 
low which antiferromagnetic ordering takes place exhibits 
a non-monotonic behavior with the ratio Uf/Wf. This 
is shown in Fig. |l^. At weak coupling T/y grows with 
Uf/Wf in a manner well-described by the Stoner crite- 
rion, and manifested in the Hartree-Fock result shown by 
the solid curve. However, at strong coupling the single 
band Hubbard model maps onto the quantum spin-1/2 
antiferromagnetic Heisenberg model, with exchange con- 
stant J — 4:t'jj:/Uf. Thus at strong coupling the Hubbard 
model Tn turns over and falls, as seen in the QMC results 
(data points) in Fig. ^ |8^]. Indeed, high temperature 
series calculations give Tn — S.SSt'jj/Uf [Q, shown by 
the dashed curve in Fig. |ri|. 

It is important to realize that Hartree-Fock (HF) mean 
field theory completely misidentifies the strong coupling 
behavior of Tn, and, instead of the correct Heisenberg 
result, predicts that Tn approaches C///4 with increas- 
ing Uf/Wf. At the HF Neel temperature, local moments 



do become well-formed, but contrary to the HF results, 
these moments remain unordered until much lower tem- 
peratures set by the scale J. The QMC calculations, on 
the other hand, can separately resolve these two different 
energy scales, either by studying the appropriate correla- 
tion functions, the local moment and the magnetic struc- 
ture factor | |76[ , or else directly from the thermodynamics 

US- 

This single band model already includes a number 
of features of possible interest to the volume collapse, 
particularly to Johansson's Mott transition picture [ p8[ . 
QMC simulations find the evolution of the density of 
electrons with chemical potential shows flat plateaus in- 
dicative of a well-formed gap as the occupation per site 
passes through (n) = 1 [Q. Unlike mean-field treat- 
ments, however, the transfer of spectral weight as the 
ratio of Uf/Wf increases is much less abrupt, and in- 
cludes the development of a resonance at the Fermi sur- 
face before the formation of upper and lower Hubbard 
bands, as sketched in the inserts in Fig. 11 |p3]. Calcu- 



lations for the infinite dimensional Hubbard model sug- 
gest a continuous metal-insulator (MI) transition in the 
vicinity of Uf/Wf ^ 1 (shaded region labeled "MI") as- 
sociated with the loss of this resonance [^Sj. The solid 
line below is an artificial extension of this boundary into 
the antiferromagnetic region, which might apply were the 
magnetic order suppressed, for example, by frustrating 
interactions. Similar infinite dimensional calculations for 
the antiferromagnetic-paramagnetic boundary are in rea- 
sonable agreement with the three-dimensional QMC re- 
sults in Fig. |ll], suggesting the relevance of this work to 
the 3~d case of interest here. 

The single band Hubbard model can, of course, be gen- 
eralized by including longer range hopping or Coulomb 
interactions. Including multiple orbitals, however, allows 
one to describe a number of phenomena, like the screen- 
ing off-electron spins by conduction bands, which occur 
in rare earth systems like Ce. The most commonly con- 
sidered two band Hamiltonian is the Periodic Anderson 
Model (PAM). 



ka k(T 

i 

+ ^efntfa - H^iriifa + nida) ■ (17) 

For the simple cubic structure considered here, we have 
taken 

efe = —2tdd [cos kxa + cos kyU + cos kzo] , (18) 
Vk = —"^tfd [cosk^a -\- cosfcyfl -|- coskza] , (19) 

where a is the lattice constant. We now have two sets of 
orbitals, d and /. The first is "itinerant" - the valence or d 
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orbitals are hybridized on neighboring sites, giving rise to 
a dispersive band efc, with no interactions, Ud = 0. The 
second is "localized," with a flat, non-dispersive band 
£/, and also highly correlated, electrons of spin up and 
down repel with an on-site energy U / . The / orbitals 
hybridize with the itinerant d band with matrix element 
Vk- There are different choices of Vk used in the litera- 
ture, the most common being a fc-independent constant 
for which the localized and itinerant electrons hybridize 
on the same site. Our choice. Eg. (|l9|) , corresponds to 
the localized orbitals hybridizing with the near-neighbor 
itinerant orbitals. 

The case e/ = is termed the "symmetric limit" of 
this model and, as in the single band Hubbard model, at 
chemical potential ^ = both the localized and itinerant 
orbitals are precisely half-filled for any temperature and 
choice of the parameters tdd,tfd, and Uf for a bipartite 
lattice. Again, there is no sign problem in this limit. Be- 
cause of the particle-hole symmetric form in which the 
third term of Eq.(|l^) is written, this limit corresponds 
to the site energy choice Sf = ef — Uf/2 = —Uf/2 in the 
notation of Fig. |^. It is our intention to eventually con- 
sider both non-symmetric limits as well as to introduce a 
dispersive / band. Nonetheless, the present simpler case 
is of some interest, especially since Fig. ^ suggests the 
/-valence interaction may be the more important origin 
of the overall /-band width. From this perspective, the 
1-band Hubbard Hamiltonian and the present Periodic 
Anderson Hamiltonian may be viewed as the simplest 
nontrivial correlated models within which to examine the 
effects of /-/ and /-valence interactions, respectively. 

The physics of the PAM is by no means entirely sorted 
out. Again, it is at half-filling where our understand- 
ing is most complete. There, a competition between an 
"RKKY" energy scale -Erkky J'^/Wd and a Kondo 
energy scale Ek oc Wdexp{—Wd/ J), where J cx fj^/Uj, 
that determines which of two possible ground states oc- 
curs. If £'rkky is dominant, at low temperatures the 
moments on the localized orbitals organize in an anti- 
ferromagnetic pattern via a coupling mediated by the 
itinerant band. If i?K is dominant, at low temperatures 
the moments on the localized orbitals form singlets with 
electrons in the conduction band. In this phase the low 
temperature Curie contribution to the uniform suscep- 
tibility is suppressed, even though the localized orbitals 
remain singly occupied. 

The phase diagram in the T/Uf-Uf /Wf plane for 
tdd = 1 eV and Uf — 6 eV is rather reminiscent of the 
corresponding phase diagram of the single band Hubbard 
Hamiltonian, where in analogy to that case we have taken 
Wf = 12tfd for the PAM, as discussed at the end of Sec- 
tion 3. As shown in Fig. there is a low temperature 
antiferromagnetic phase, whose transition temperature 
exhibits a non-monotonic dependence on Uf/Wf similar 
to the one band case. There is also a cross-over temper- 
ature which separates the paramagnetic phase into a 




U/W, 

FIG. 12. Phase diagram of the 3-d, symmetric Periodic 
Anderson Model (PAM) at half-filling (1 electron per or- 
bital per site) from QMC simulations. Here Uf = 6tdd 
is fixed, Wf — 12tfd is varied, and the number of sites 
7V = 4x4x4 = 64. QMC data points are shown by tri- 
angles, the dark triangles indicate antiferromagnetic order- 
ing in the f-band. The antiferromagnetic region is indicated 
by shading. At higher temperatures and high Uf/Wf free 
f-band moments still exist, but are decoupled and not or- 
dered. As Uf/Wf decreases, the f-band electrons begin to 
form singlets with the conduction-band electrons, until in the 
regime labeled "Kondo singlets" the Kondo-like singlets are 
fully-formed. The hashed region marks the Kondo singlet 
formation cross-over. Experimentally the two regimes are 
distinguished by the absence of a Curie-like 1/T divergence 
of the uniform susceptibility in the Kondo regime. Hartree 
Fock and RKKY curves are shown for comparison. 

region with decoupled moments (to the right of the 
hatched boundary), and a region where these moments 
form singlets with electrons in the conduction band (to 
the left). This cross-over region is the counterpart of 
the metal-insulator transition in Fig. in that in both 
cases this boundary is associated with appearance of the 
central resonance for decreasing Uf/Wf. The Hartree- 
Fock prediction for Tn (solid curve) in Fig. |l^ also has 
the same appearance as in Fig. first rising and then 
saturating at Uf/4 (not shown) for increasing Uf/Wf. 
The RKKY energy, EnKKv/Uf cx (Uf/Wf)-^ for con- 
stant Wd and Uf, is analogous to the Heisenberg curve 
in Fig. |ll|. The ratio Srkky/C^/ is plotted as the dashed 
curve in Fig. |l^, and has been scaled to be consistent 
with the QMC location of Tn near Uf/Wf = 0.6. 

The parameters in Figs, ^-j^ may be used to translate 
the location of the observed rare earth volume collapse 
transitions (or symmetry changes in Nd and Sm) into val- 
ues of Uf/Wf. Considering only the /-valence hopping, 
these are in the range Uf/Wf = 1.1-1.9 except for Ce 
which is 3.7-2.6, as noted earlier. It is not unreasonable 
to associate these values with the cross-over region just 
discussed, Uf/Wf ^ 0.8 or larger if extrapolated to room 
temperature, whereas by contrast, the Hartree-Fock pre- 
dictions for the symmetric two-band model indicate a 
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transition for Uf/Wf ~ 0.2. It should also be empha- 
sized that the one and two-band models discussed here 
correspond to half-filled / bands, and that lower /-band 
filling should favor the "itinerant" states, pushing tran- 
sitions to larger values of Uf/Wf, as may occur with Ce. 

C. Density of States in the PAM 

The behavior of the quasiparticle density of states 
(DOS) N{uj/Uf) provides an especially graphic illustra- 
tion of how QMC differs from mean field theories such as 
HF in its treatment of correlations. 

In Fig. |l^ we show Nioj/Uf) of the localized band 
of the PAM for different inverse hybridizations Uf/Wf 
at a temperature above the Neel temperature Tjv at 
which antiferromagnetic ordering occurs. The DOS 
was obtained by numerical analytic continuation of the 
single-particle imaginary-time Green's function, G{t) = 
W I]p('^rCp(T)cj,(0)}, using the relation 

poo -TUj/Uf 

G{t) = / d{Lu/Uf) N{Lu/Uf) (20) 
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FIG. 13. The f-band quasi-particle density of states N{ij) 
as a function of distance from the Fermi energy {uip = 0) 
of the PAM for T/Uf = 1/30, Uf = Gtaa, and at half-filling 
(n) = 2. Two values of Uf/Wf are shown. For Uf/Wf large, 
the two Hubbard-band peaks contain all of the weight in the 
DOS (solid line). [The solid line values have been reduced by 
a factor of 3.] As Uf/Wf decreases, the Kondo-singlet forma- 
tion regime is entered and spectral weight shifts continually 
to the center, Fermi-energy region. The dot-dashed line rep- 
resents a small Uf/Wf where the singlets are well-formed. 
The pseudo-gap, or reduction in the spectral function N{uj) 
at the Fermi energy, is due to the insulating nature of the 
symmetric, half-filled PAM at small Uf/Wf. 



for imaginary time t, temperature T/Uf, and frequency 
uj/Uf. We use the maximum entropy method to 
perform the analytic continuation in a technique which 
utilizes the full imaginary-time covariance matrix. 

For a large Uf/Wf (sohd line), the spectral weight is 
concentrated in the lower and upper 'Hubbard' bands, 
corresponding to singly-occupied localized orbitals and 
to doubly-occupied localized orbitals, respectively. As 
Uf /Wf is decreased, the hybridization between the bands 
increases and the spectral weight simultaneously begins 
to shift continuously from the two peaks to the region 
near the Fermi energy {uj/Uf = 0). The build-up of the 
central peak in the DOS corresponds to the formation of 
singlets between electrons on the localized and itinerant 
orbitals as will be discussed shortly. At a sufficiently 
small Uf/Wf a fully-developed three-peak structure is 
seen in the DOS (dot-dashed line). The pseudo-gap, 
or reduction in the spectral function N{uj) at the Fermi 
energy, is due to the insulating nature of the symmetric, 
half-filled PAM at small Uf/W. 

There is no mean-field analogy either for the three 
peak structure in the DOS, or for this continuous shift 
in the spectral weight. At these same conditions, HF 
yields a partial /-DOS similar to Fig. ^for the localized 
limit, Uf/Wf = 5. However, by Uf/Wf = 0.5, these 
same lower and upper Hubbard peaks are still present, 
although considerably broadened. Not mitil Uf/Wf = 
0.22, at T/Uf = 0.033, do the HF antiferromagnetic and 
paramagnetic free energies cross, at which point a single 
broad central HF /-DOS peak is observed for smaller 
Uf/Wf. 

D. Magnetic Properties of the Periodic Anderson 
Hamiltonian 

The ultimate goal of using QMC to compute precisely 
the contribution of the correlation energy for the volume 
collapse transition is still in the future, both in deter- 
mining the specific appropriate lattice Hamiltonian and 
in then carrying out the computations. Here we will show 
some further results which illustrate that QMC can sen- 
sitively pick up the basic physics of correlations between 
moments like Kondo singlet formation and long range an- 
tiferromagnetic order. We will continue to simulate the 
PAM, Eq.(|l|) for tdd = 1 eV, Uf = 6 eV, and a range 
of tfd values indicated as before by Wf = 12 tfd. These 
results will provide some of the background to the phase 
diagram shown in Fig. |l^. 

While we are primarily interested in determining the 
thermodynamics to infer the equation of state, it is useful 
to look at magnetic properties such as the local moment 
and magnetic structure factor. In Fig. |l^ we show the 
square of the local moment A'd(/) = {^id(f) ■ ^id(f)) on the 
itinerant {d) and localized (/) orbitals. ft^ takes on a 
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FIG. 14. Square of the local moment on the itinerant (lower 
curve) and localized (upper curve) orbitals at various temper- 
atures T/Uf and Uf/Wf. The values shown have the total 
moment for a non-interacting system subtracted, so that on 
the curves corresponds to the moment of a non-interacting 
electron system. The itinerant orbital electrons remain at 
approximately the non-interacting moment value. The local 
(f) orbitals have increasingly well-formed moments as Uf/Wf 
becomes larger. 

value essentially equal to that for a non-interacting elec- 
tron gas. This is a consequence of the lack of correlations, 
Ud = 0, between up and down spin electrons. In contrast, 
/Xj shows a significant enhancement throughout the range 
of interband hybridization tfd. 

jjL^j is the most local measure possible of magnetism- 
the presence or absence of a moment on a site. The 
magnetic structure factor determines the degree of the 
order between spins on different sites across the entire 
lattice, 



(21) 



This quantity is shown in Fig. ^ for various values of 
the temperature and interorbital hybridization, tfd^ re- 
flected in the ratio Uf/Wf. A careful finite size scaling 
study is necessary before rigorous conclusions concern- 
ing the presence of long range antiferromagnetic order 
are drawn. The above data, however, suggest a strong 
tendency for the moments to order antiferromagnetically, 
with an ordering temperature which is maximal around 
Uf/Wf « 0.6. 

As discussed above, singlet formation between the 
spins of the localized and itinerant electrons competes 
with antiferromagnetic order. Fig. 16 shows a quantity 
which measures this tendency. Here, Cfd is essentially 
(Si/ • Sjd-) averaged over near-neighbor pairs, i, j, and 
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FIG. 15. Anti-ferromagnetic structure fac- 

tor, S'//(r/(7/,q* — (7r,7r,7r)) for various values of Uf/Wf. 
The horizontal line demarks the existence of antiferromag- 
netic long-range order in the system. 

appropriately normalized to reflect the number of such 
pairs and the average on-site local moments. The cross- 
ing of these curves through 40-60 percent of their lowest 
saturated value, indicated by the horizontal lines, gives 
an estimate of the Kondo temperature below which sin- 
glet formation takes place [^ . 

From the results shown in Figs. |r^[l6|, and analysis of 
other quantities such as the density of states, wc infer the 
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FIG. 16. A measure of the tendency for Kondo singlet for- 
mation. The horizontal lines indicate the approximate range 
of the cross-over between the paramagnetic and Kondo sin- 
glet phases (hatched region of Fig. [l^ .) 
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magnetic phase diagram in the {Uf/Wf, T/Uf) plane for 
the PAM on a 3-d simple cubic lattice at tdd = 1 eV, 
Uf — 6 eV, and half-filling, shown in Fig. p^ . 

Low-temperature HF calculations of the squared lo- 
cal moments, A'^(j), resemble the QMC results in Fig. |l^ 
for the antiferromagnetic phase. However, neither f or 
d moment shows any enhancement in the paramagnetic 
phase, demonstrating the well known fact that moment 
formation and magnetic order coincide in HF. Both anti- 
ferromagnetic and paramagnetic solutions exhibit values 
of (Si/ ■ S2d) which decrease from ~ to ~ —1/4 as 
Uf/Wf is reduced, where 1 and 2 are the two sites of 
a doubled unit cell. While a crude immitation of the 
Kondo-like behavior in Fig. ^ even at T = 0, this varia- 
tion is quite gradual, and the mean-field expectation can 
of course never approach the full singlet limit of —3/4. 

E. Thermodynamic Properties of the Periodic 
Anderson Hamiltonian 

The energy can be easily measured as a function of 
temperature in QMC. From this data, the specific heat 
can be obtained, for example using a finite difference ap- 
proximation for the derivative at two relatively closely 
spaced temperature values. A more useful way to proceed 
is to measure E{T) on a sufficiently fine grid of temper- 
atures and fit to a physically motivated functional form. 
We have chosen 

N 

E{T) = E{0) + c„e-''"^ , (22) 

n=l 

where A and c„ are free parameters. N is chosen to be 5- 
6 to compromise between allowing sufficient flexibility to 
pick up the full structure without ovcrfitting. A check on 
this procedure is provided by computing the integrated 
area, 

^~,r^fl = f;i|^41„2-S.. (23) 

■^^> n=l 

Here the presence of the term 5*0 is a result of the use 
of a finite periodic cluster where the entropy does not 
go to zero at T = |^^. We find this sum rule on 
the total entropy is obeyed to within less than 2 percent 
for tfd > 0.8. For smaller tfd the system orders anti- 
ferromagnetically at temperatures below those accessed 
in the simulations (T > 0.1 eV), and the integral instead 
approaches 3 In 2 — iSq , correctly reflecting the missing en- 
tropy of spin ordering. Note that the QMC calculations 
were extended to sufficiently high temperatures to estab- 
lish the large T asymptotic dependence of C{T)/T. We 
have also tried other fitting forms, and find our results for 
the specific heat are rather insensitive to the particular 
functions involved. 




FIG. 17. Specific heat of the PAM obtained by fitting the 
QMC energies. The area under the curves satisfy the sum 
rule, Eq. bs. The peaks in the specific heat appear to be as- 
sociated with antiferromagnetic ordering. (See, also, Fig. ^) 
The corresponding peaks at tfd = 0.2 and 0.4 lie below 
T < 0.1 eV, the temperature range sampled by QMC, and 
so are absent. 

The results of this procedure are shown in Fig. |l^. 
The specific heat exhibits a broad bump at high temper- 
atures, the energy scale Uf of charge fluctuations, that 
is, excitations to the manifold of doubly occupied states. 
At low temperatures there is a much sharper peak which, 
given the data on the magnetic correlations, we believe 
is associated with antiferromagnetic ordering. We do not 
see any sharp feature associated with singlet formation, 
which drives the volume collapse transition according to 
the Kondo picture ||3^,^. However, it is well known 
that the signature of singlet formation in the specific 
heat becomes sharp only for relatively small tfd, where, 
in this symmetric half-filled model, antiferromagnetism 
strongly competes 

The resolution of this difficulty is to study the PAM 
in parameter regimes where the antiferromagnetic order 
is suppressed, either by going away from the symmetric 
limit ef = or else away from half-filling /i = 0. As 
discussed above, doing so is also suggested by the LDF 
analysis, where the / site energy, Sf = e/ — C///2, is 
seen in Fig. to rise with compression relative to the 
valence chemical potential. Moreover, HF calculations 
show there is a significant suppression of the Neel tem- 
perature by changing e/ or fx. We have found this result 
using the difference in HF energies between the antifer- 
romagnetic and ferromagnetic solutions as an estimate of 
Tn. This approximation gives the correct strong coupling 
Heisenberg scaling T/v « t'^ /U in the single band model, 
as opposed to T/v oc U which is the formally correct HF 
prediction obtained from looking at the temperature at 
which the free energy of the antiferromagnetic solution 
becomes lower than the paramagnetic one. 
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HF calculations of C{T)/T for the same 4^-site two- 
band model bear only a gross resemblance to Fig. |l^. 
There is the appearance of a magnetic ordering peak 
caused by the drop in C (T) /T from antiferromagnetic to 
paramagnetic solutions with increasting T at T/y. How- 
ever, this occurs at the HF T^r which incorrectly ap- 
proaches J7//4 = 1.5 eV, rather than 0, as tfd — > 0. As- 
suming the same band widths, Wfd, the rare earth transi- 
tions would occur for values of tfd in the range 0.1-0.3 eV 
for near-neighbor interactions in the simple cubic struc- 
ture. At tfd = 0.2 eV, the HF paramagnetic C{T)/T is 
in reasonable agreement with the corresponding curve in 
Fig. |l^ except for having a larger area (4 In 2—25*0 versus 
3 In 2— 5*0). In the experimentally observed paramagnetic 
regime, the HF paramagnetic calculation provides a bet- 
ter approximation to the true thermodynamics than does 
the HF antiferromangnetic solution even though the lat- 
ter may have lower free energy. 

VI. SUMMARY 

This paper has discussed our current experimental and 
theoretical understanding of the volume collapse transi- 
tions in the trivalent rare earth metals. The data ex- 
hibit "localized" 4/-electron behavior at low pressure, 
characterized by isolated-ion-like magnetic moments and 
the absence of these electrons in the crystal bonding, 
and "itinerant" 4/-electron behavior at higher pressures, 
without apparent moments but with clear evidence of 4/ 
bonding participation. The two regimes are separated by 
first order phase transitions which in three cases (Ce, Pr, 
and Gd) are accompanied by unusually large (9%-15%) 
volume changes. It is generally believed that these vol- 
ume collapse transitions are associated with change in the 
degree of 4/ electron correlation. Consistent with this 
idea, we have reported calculations of relevant param- 
eters which indicate the transitions occur in the range 
Uf /Wf — 1-2, where C// is the onsite Coulomb repulsion 
between 4/ electrons, and Wf is a measure of the 4/ 
band width including both /-/ and /-valence hybridiza- 
tion effects. The one exception is Ce, where the tran- 
sition occurs at a larger (^^ 3) value of this ratio which 
may reflect the low (/^) /-band filling in this case. 

We have compared a variety of orbitally realistic, 
mean-fleld based methods used to calculate these transi- 
tions: spin and orbital polarized local density functional 
theory, self-interaction corrections, the LDA+U method, 
and Hartree-Fock (HF) with static screening. As a class, 
these methods represent the localized phases by polar- 
ized solutions with an integer number of occupied, split- 
off 4/ bands, and the itinerant phases, by a single broad 
collection of 4/ bands overlapping the Fermi level. The 
transition is therefore accompanied by an abrupt trans- 
fer of spectral weight (the split-off bands) to the vicin- 
ity of the Fermi level. Both the loss of moment and 



the onset of 4/ bonding are clear from these calcula- 
tions. At very small and large volumes, Uf « Wf 
and Uf >> Wf, respectively, such mean-field solutions 
should provide the correct ground states, aside from mis- 
treatment of intraatomic correlation for multi-/ electron 
atoms in the latter limit. Indeed, Anderson has used 
the HF solution to elucidate the limits of the one-band 
Hubbard model [Q. The problem hes in between these 
limits where competition between the two energy scales 
leads in general to a correlated solution. Thus, while 
several calculations have reported reasonable agreement 
with experiment for the location of the transitions ||5^]- 
|6l| , there are substantial differences between HF and ex- 
act Quantum Monte Carlo (QMC) solutions for few band 
Hamiltonians which raise doubts about these approaches. 

To assess correlation effects, we have reviewed exact 
QMC solutions for the one band Hubbard model, and re- 
ported new QMC results for the two-band periodic An- 
derson lattice model. We contrast both solutions with 
HF results for the same Hamitonians. The QMC results 
for both one and two-band models show non-monotonic 
behavior of the Neel temperature, T/v, as a function of 
Uf/Wf, where this ratio scales like volume to a power 
^ 2. T/v first rises with increasing Uf/Wf as antifer- 
romagnetic order accompanies moment formation, and 
then falls as decreasing exchange interactions lead to mo- 
ment disorder. HF misses the latter effect, so that its T/v 
continues to rise to a large saturation value of J7//4 ~ 1.5 
eV. The HF transitions are therefore always intrinsically 
associated with both moment formation and magnetic 
ordering. The region of magnetic order is greatly sup- 
pressed in the QMC solutions, e.g., lying below 0.1 
eV for the present PAM results based on realistic param- 
eters. This value can be further suppressed by chang- 
ing the / level position (non-symmetric PAM) , consistent 
with the experimental observation that the volume col- 
lapse transitions occur entirely within the paramagnetic 
regime. 

The dramatic new feature in the correlated QMC so- 
lutions, which has no counterpart in mean field, is the 
three peak structure in the / density of states (DOS) 
found especially at temperatures just above the region of 
magnetic order. The growth of the central Fermi-level 
resonance for decreasing Uf /Wf comes from a continu- 
ous transfer of spectral weight from the outer split peaks, 
and is accompanied in the QMC results with reduction 
of the magnetic moment as might be deduced from the 
magnetic susceptibility. This reduction occurs at consid- 
erably larger values of Uf/Wf or volume than docs the 
HF transition, and is in much better agreement with our 
"experimental" Uf/Wf placement of the collapse transi- 
tions. 

The central resonance is the well known signature of 
the Kondo effect for the PAM, and in its impurity An- 
derson model context, the critical aspect of the Kondo 
volume collapse model for Ce p9|,Eo|. However, it is in- 
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triguing that the infinite dimensional one-band Hubbard 
model reduces formally to the impurity Anderson model 
and exhibits a similar central resonance, which might be 
associated with / electrons screening their own moments 
p3[ . This blurs the distinction between the Mott transi- 
tion picture [Q, traditionally associated with the Hub- 
bard model, and the Kondo volume collapse |3^,Q . The 
distinction is further complicated in any mean field ap- 
proach which would pick up only half (e.g., fid^ but not 
of a Kondo like singlet, and might therefore still 
look like / magnetic order. This is not the case for QMC 
calculations which can separately identify reduction in 
the local moments as well as their screening by the va- 
lence electrons, as seen in the present PAM results. More 
definitive use of QMC in this manner, however, awaits 
studies in other parameter regimes (e.g., non symmetric 
PAM) in search of thermodynamic features, for example, 
which might unambiguously demonstrate that the given 
model does indeed capture the essential nature of the rare 
earth volume collapse transitions. 

A natural future direction for the theory of the rare 
earth volume collapse transitions would be to compro- 
mise between the all-orbital mean-field based methods 
and few-orbital exact treatments, so as to incorporate 
only the most important correlation effects within a suf- 
ficiently rich multi-orbital framework. There appears to 
be, for example, a growing awareness that some form of 
dynamic screening is critical to incorporate features such 
as the three peak density of states, requiring inclusion 
of an energy-dependent self-energy. Some success has 
been achieved with local, k-independent approximations 
to such self energies [89 -91 . In principal this approach 



reduces to a correlated impurity problem which could be 
solved by QMC methods for far more orbitals per site 
that possible in the fully correlated lattice. Should in- 
termediate or longer-range correlations be critical, a new 
modified dynamical mean field technique may allow in- 
clusion of such correlations in addition to multiple bands 
p2| . All of these techniques form a bridge between local 
density functional methods which serve to define realistic 
multi-orbital Hamiltonians which may then be solved by 
correlated approaches such as QMC. 
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